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SPACE SHUHLE POGO ACTIVE CONTROLLfR DESIGN USING 
FREQUENCY DOMAIN OPTIMIZATION 
by Robert C. Seidel, Carl F. Lorenzo, and Bruce Lehtinen 
Lewis Research Center 

SUMMARY 

A frequency domain parameter optimization technique was used in a preliminary de- 
sign study of active pogo suppression controls for the space shuttle. The technique uses 
a conjugate gradient search procedure to optimize a cost fimction of controller perform- 
ance and is well suited for designing low -order controllers for higher order systems. 

The shuttle model was a two -pump and six-structural -mode linear model of order 41, 
representing a worst-case condition. The suction pressure of the high-pressure oxidizer 
pump was the controlled variable. The parameter optimization was performed for four 
different low-order controllers. The resulting closed -loop system performances were 
evaluated by several means including controller performance averaged square values, 
frequency responses, system damping ratios, and sensitivities to structural model 
changes. A promising feedback controller structure was found to be a lead -lag design. 


INTRODUCTION 

Liquid -propellant rocket vehicles can experience oscillations of thrust commonly 
known as pogo. The pogo problem results from interactions between the vehicle struc- 
tural modes of vibration and the propulsion system modes. Two concepts applicable to 
the space shuttle suppress pogo by suppressing feedline flow changes. They are a pas- 
sive accumulator and an active suppressor. The active suppressor consists of a piston 
type of device which injects and removes flow from the propulsion system based on a 
sensed signal (ref. 1). A critical element in an active suppressor is the feedback con- 
troller design which is the subject of this report. 

The controller design uses a frequency -domain parameter optimization technique 
that is well suited for designing low -order controllers for higher order systems (ref, 2). 
The pogo model of the space shuttle is basically the Aerospace Corporation linear model 



of references 3 and 4, At Lewis the model was cast into a matrix state variable model 
with multiple structural modes and lumped feedline description. Plant transfer functions 
were obtained from the matrix model using the computer programs of reference 5. 

This report describes the control problem and discusses the design results. The 
control problem description includes a discussion of the shuttle model, the controller, 
and the parameter optimization technique. The results present controller performances 
(averaged square values) for four controller structures. Analysis of selected designs 
includes system frequency responses, eigenvalue damping ratios, and sensitivities to 
changes in structural mode natural frequencies. 


SYMBOLS 


b controller parameter vector 

bj b vector element proportional gain, dimensionless 

\>2 b vector element derivative gain, dimensionless 

bg b vector element lead -lag gain, dimensionless 

b^ b vector element lead -lag zero, rad/sec 

bg b vector element second-order lead-lag gain, dimensionless 

bg b vector element second-order lead -lag damping ratio, dimensionless 

bry b vector element second-order lead-lag natural .frequency, rad/sec 

C cost function, (N/cm^)^ 

D power spectral density of d, [cm /(sec* rad)] 

Q 

d disturbance in feedline, cm /sec 

2 2 

F cost function value, (N/cm ) /(rad/sec) 

2 3 

G system open-loop transfer function, (N/cm )/(cm /sec) 

g 

actuator, cm /sec 

Gel system closed -loop transfer function, (N/cm )/(cm /sec) 

H controller, (cm^/sec)/(N/cm^) 

h controller filter 

jo; Fovu*ier variable, 1/sec 

K pressure transducer gain, cm /N 


2 


2 

actuator gain, cm /sec 

2 

P,^ HPOP suction pressure variation, N/cm 

Q 

Q„ suppressor flow, cm /sec 

a 

2 3 2 

r weighting factor [(N/cm )/(cm /sec)] 

s Laplace variable, 1/sec 

t time, sec 

2 2 

V power spectral density of v, [N/(cm -rad)] 

2 

V measvirement noise, N/cm 

y variable 

T random noise averaging time, sec 

o) angular frequency, rad/sec 

suppressor actuator pole, rad/sec 

A 

Subscripts: 

d excitation by d 

( normalized value 

V excitation by v 

Superscripts: 

_ frequency average 

- time average 


CONTROL PROBLEM DESCRIPTION 
Description of Shuttle Model 

A drawing of the space shuttle ca!rrying the orbiter vehicle is shown in figure 1(a). 
The liquid- oxygen feed system dynamics can interact with the vehical structural dynam- 
ics to produce pogo. In figure 1(b) a schematic diagram is shown of the liquid oxygen 
part of the shuttle propulsion system that includes the liquid oxygen tank, feedline, low- 
pressure oxidizer pump (LPOP), the suppressor (under consideration herein), high- 
pressure oxidizer pump (HPOP), and the rocket chamber. 

To obtain the (state variable matrix) model used in this report, two changes were 
made to the reference 3 model. The first change was to approximate the distributed 
parameter feedline equations with a number of lumped parameter approximations. The 
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feedline segment between stations 1 and 2 was divided into nine parts, and the segment 
between stations 2 and 3 was divided into three parts. The second change was to include 
six structural modes in the model simviltaneously instead of only one mode, as was done 
in reference 3. This insured that a controller designed to stabilize one mode had no se- 
rious destabilizing effects on the remaining modes. A detailed discussion of the mod- 
ified model is presented in the appendix. 

The model of reference 3 has variable model parameters that change for different 
flight conditions. A nominally worst-case model was assumed for this control study to 
avoid studying all possible configurations in flight. The worst -case model is basically 
the reference 3 model at the end -burn conditions except for the following: the end -burn 
structural mode at 2.81 hertz was replaced by a less stable configuration which was the 
2. 3 -hertz mode from the after -separation condition. The nominal 26-hertz structural 
mode was shifted to 24 hertz to make the system open-loop unstable. Also, the low- and 
high-pressure oxidizer pump gains m^ + 1 and mg + 1 were set to the maximum ex- 
pected value (over the flight time) of 2.20 and 1. 54, respectively. The pump inlet bubble 
compliances were set to the minimum (end burn) values of = 5.944 centimeter per 
newton (0. 25 in. ^/Ibf) and = 2. 376 cm^/N (0. 10 in. ^/Ibf). 


Description of Control System 

The control system structure assumed as a framework for the controller design 
problem is shown in figure 2. The plant is block G(s) (where s is the Laplace variable) 
and the plant output is P,^, a pressxire perturbation from the nominal steady -state oper- 
ating point of the HPOP suction pressure. This pressure variation is sensed and used 
as the system feedback signal. The suppressor considered in this study was located at 
the HPOP inlet. The proximity of the pressiire measurement to the suppressor reduces 
plant phase lag and will generally resvdt in better closed-loop control. The plant dis- 
turbance d is assumed to be a flow disturbance that enters the plant in parallel with the 
system suppressor flow (fig. 2) at station 7. This is compatible with the optimiza- 
tion program constraints. There is also a measurement noise signal v. The d and v 
signals are assumed to be independent Gaussian white noises with power spectral den- 
sities D and V, respectively. 

The controller H(s, b) representation consists of three factors; a transducer gain 
K, a suppressor actuator transfer function G (s), and a controller filter transfer func- 

A 

tion h(s, b). The filter includes a vector of variable parameters b to be determined as 
described in the section Parameter Optimization. 

H(s,b) =Kh(s,b)G^(s) (1) 
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The suppressor actuator transfer function G.(s) was assumed to be a first-order lag 

ci 

with a pole u} and a steady -state gain K . 

3 . cL 


G^(s) 


K. 


- 1+1 




( 2 ) 


For the purposes of this study the following values were chosen for the controller 
constants: 

K = 1. 450 cm^/N (1. 000 in.^/lbf) 

K = 16. 39 cm^/sec (1. 000 in. ^/sec) 

3 

= (27t) 50 rad/sec 

3 

A more detailed design study should consider an actuator representation including higher 
order servovalve characteristics and saturation limits. 

The purpose of the controller is to minimize the variation in the system pressures 
and flows due to the disturbance d. At the same time it is desirable to minimize the 
control effort to accomplish this Q^. The pressure variation P,^ has been used in this 
study as an indicator of the system response. 

The plant -transfer -function G(s) = P,^(s)/d (of fig. 2) frequency -response magnitude 
and phase is plotted in figure 3. The response is for the multistructural-mode, worst- 
case model previously described. There are a number of resonances indicated by the 
magnitude and phase data. These include a very lightly damped mode at about 2. 3 hertz 
and an unstable mode at about 24 hertz. The instability is not apparent from the separate 
frequency response plots, but the system can be seen to be unstable from an eigenvalue 
analysis. 

The complex roots of the denominator of the plant transfer function set the plant 
resonances. The damping ratios and natural frequencies associated with the complex 
roots (the eigenvalues of the system matrix) are tabulated in table I. The natural fre- 
quencies associated with the uncoupled structural mode natural frequencies are footnoted. 
The uncoupled structural modes are assmned to have a damping ratio of 0.01 in the 
model. The coupling of the structiiral and prop\ilsion system modes has reduced some 
system eigenvalue damping ratios. The mode at about 2.3 hertz has a damping ratio of 
only 0. 0008. The unstable mode at about 24 hertz has a damping ratio of -0. 0015. 

The frequency range of major interest for the study is below 30 hertz as set by the 
highest structural mode at 26. 8 hertz. There are low damping ratios for the frequencies 
above about 40 hertz. But these higher frequency modes are associated with the feedline 
and do not involve a pogo problem with the structural modes. The approximations made 
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to form a lumped parameter model of the feedline dynamics lose accuracy as frequency 
is increased. 


Parameter Optimization 

A cost function, defined here, reflects the average variations in P,y and 
caused by the random disturbance and measurement noises d and v (fig. 2). The cost 
function C formulated was the mean square value of the pressure variation of P,^ plus 
the mean square value of the control effort weighted by scalar r. 

/ T = ■ ■■ 

[P,(t,b)2 + rQ^(t,b)2]dt = P^(b)2 + rQ^(b)2 (3) 

-T 

where the double bar notation denotes a time average. By Parseval’s theorem this cost 
function can be transferred out of the time domain and into the frequency domain where 
the parameter optimization calculations are made: 

C(b) = i f [|p^ac^,b)|^+r|Q^0a),b)|2]da> = l |P^(b)|2 + r|0^(b)|2] (4) 

Jq it 

where joj is the Fourier variable, w is angular frequency, and the bar notation denotes 
integration with respect to frequency. Because the inputs are independent and the system 
is linear, equation (4) can be expanded to express the individual effects of disturbance d 
and measurement noise v as, 

C(b)=i |P,(b)||* |P,(b)|2 + r|Q^(b)|2+r|Q^(b)|2 (5) 

7T L 

where the subscripts d and v indicate the origin of the excitation. 

As was done in reference 2, C is expressed in terms of transfer function relation- 
ships. The expression for C is obtained in terms of the following known quantities: the 
power spectral density D of the disturbance d, the power spectral density V of the 
measiu'ement noise v, the plant transfer function G, and the controller transfer function 
H. To start, the system closed -loop transfer function Gel between P,^ and d is the 
usual closed-loop feedback relationship: 
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( 6 ) 


Ggj(jw,b) = 


From block diagram algebra, we obtain 


G(jg)) 

1 + G(iw)H(ia»,bl 


|p,^(ja,,b)|2.D|G^l(ja),b)|2 

I ^ 7(3 w,b ) \l = V| G^j(jo),b)- H(jo),b) 1 2 

iQaO^^.b)!! = D|Gci(i‘^,b)*H(ja),b)|2 






(7) 




= V|G^j(jo),b)-H(ja),b)/G(jco)| 


and C(b) may now be expressed as 


C(b) 


---f 

TT -^0 


F(ct»,b)do) 


where, after substituting and combining terms, F is 

■wn / j \ I /* , \ (2 


f 


1 

||H(jo;,b )|2 

rU + V + 


|G(jw)|^_ 

J 


( 8 ) 


(9) 


With an equation from which to evaluate C as a function of known system relation- 
ships and variable controller parameter initial estimates b, the parameter optimization 
problem is to minimize C(b) by the proper selection of parameter vector b. Most 
parameter optimization methods, including the conjugate gradient method (which is used 
here) require also that the gradient of C(b) with respect to the parameter vector b be 
calculated. The gradient of C(b) is obtained by taking the gradient of both sides of equa- 
tion ( 8 ). 

VC(b)=i / VF(co,b)dw (10) 

n -'0 


Use is made of the identity 
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v|y(joj)|^ = 2 [y(-jw)Vy(ja))] 


and equation (9) to express VF(o),b) as 


/ 


1 

VF(o;,b) =2|G^l(jaj,b)|2/?eh 

H(-jw,b) 

rP + V + 

|G(jo,)|2_ 


- DG(ju)) 


1 G^j(j«,b) 




G(jo)) 


VH(ja),b)| 


( 11 ) 


Once the structure for H(joj,b) has been chosen and the values of b selected, the cost 
function and its gradient can be calculated. 

The computer program which computes the cost function and the cost function gra- 
dient and which conducts the gradient search is basically the program from reference 2. 
The trapezoidal integration rule is used to evaluate the cost fimction and gradient inte- 
grals (eqs. (8) and (10)). The integration requires a (single) previous evaluation of the 
plant transfer function G(jo)) at discrete frequency points. For the pogo problem 219 
points were selected over a frequency range of 0. 5 to 114 hertz. The spacing of the fre- 
quency points is determined by the change in the open-loop plant, so that integration er- 
rors are minimal. In addition, a variable step size, incorporated in the reference 2 
program, further reduces errors caused by a fixed step size when the closed-loop trans- 
fer fxmction magnitude |Ggj(j,b)| exceeds the open -loop transfer function magnitude 
|G(jco)| . The additional closed -loop points calculations are based on the interpolation of 
the neighboring open-loop data points. The controller parameters b in h(s,b) are 
handled in factored form. A scaling of each parameter by its own magnitude leads to 
rapidly converging searches. The value of the cost function does not necessarily indi- 
cate whether the system is stable for the controller foimd by the search. System sta- 
bility is determined by a Routh-Horwitz stability test on the system characteristic poly- 
nomial 1 + G(s)H(s,b) (using a subroutine obtained from ref. 6). 

The computation time for a single evaluation of the cost function and gradient for a 
two-parameter controller was about 0. 36 second for the computer used (IBM 360-67 TSS). 
A typical two-parameter search would usually converge in less than 10 iterations; thus, 
a parameter optimization could usually be performed in several seconds. 


RESULTS 

Controller Filter Structures 

The best structure for the controller transfer function h(s,b) is not known on start- 
ing a controller design. The best structure may be a compromise between low order for 
reasons of practical implementation and high order for improvement in control. The 
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controller structures reported here start at order zero and increase in complexity to 
order two. Four controller structures are considered; 

(1) A proportional gain controller for which h(s,b) = 

(2) An approximate derivative controller for which 



(3) A lead -lag controller for which 


h(s,b) 



(4) A second -order lead -lag controller for which 


h(s,b) 





The fixed poles were numerically set to 1000 radians per second to obtain a realizable 
structure. If the poles were not fixed, the search tended to drive the poles toward in- 
finity with only a minimal improvement in control. 

It is interesting to note that passive accumulators could be sized using this program. 
This is true so long as the device senses local pressures and modifies local flow and the 
defining equations are linear. The derivative controller could be considered a passive 
compliant accumulator because the defining equations are essentially the same. 


CONTROLLER DESIGN 
Performance Comparisons 

Controllers were designed for each of the four stnictures presented in the previous 
section. The designs are for the case of no measurement noise. The weighting variable 
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r was changed to achieve different control efforts and performance. The resulting con- 

2 

troller designs are compared in figure 4 on the basis of the design variables P,^ and 
0? averaged square values. Each point on a curve corresponds to a parameter opti- 
mization for a particular value of the weighting variable r. The values of the design 
variables are reported as normalized values. The normalization of the pressure error 
P,y was by the reference value obtained by integrating the square of the open-loop fre- 
quency response magnitude. The control level Q„ was normalized by using the level 

a 

obtained from the derivative control structure with b2 = 80; that is, 


and 


where 


and 



7 1 ref. 



N 


12 


7 1 ref. 


= 1.5 


^cm 


cm 


cm ’ 
i sec* rad; 


Uec* radJ 


Ibf 


= 3.2 


<in. 


m. 


Isec'rad; 


D| 


f in. ^ 
Uec* radi 


Qalref. = 200 000 D 


cm 

,sec* rad; 


= 710 D 


in. 


isec*rad 


2 
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The four controller structure curves (fig. 4) all show a decrease in pressure varia- 
tion with an increase in control level (decreased r weighting). The first- and second- 
order lead -lag controllers reduce the pressure error more than the proportional con- 
troller over the range of control effort. Except for the case of proportional control, the 
pressure error is driven to vanishingly small values with increasing control effort. The 
derivative and lead-lag controller designs develop an almost vertical slope. (The pro- 
portional design also dropped, but at a greater control effort than shown in fig. 4. ) The 
sharp drop indicates that very large controller gains reduce the pressure error to very 
small values without driving the system unstable. Use of such high gains would probably 
not be permitted if the suppressor actuator was more accurately modeled with a higher 
order transfer function and saturation limits were considered. 

The curve for the derivative control appears as a shorter and more vertical line. 
The derivative control required about 84 percent of unity normalized control effort (and 
a minimum bg gain of about 0. 65) to achieve system stability. It was the plant mode at 
about 2. 3 hertz that made the system unstable for lower derivative gains. The propor- 
tional control was able to achieve a stable system with only about a fourth of the unity 
normalized control effort. For the proportional control case the plant mode at about 
24 hertz was the limiting factor . 

The lead-lag controllers were able to achieve stability with minimum control level. 
System stability for these controllers was achieved with less than about 5 percent of 
unity normalized control effort. The reduced effort is of course reflected back into the 
hydraulic flow requirements associated with a given controller design. 


Controller Parameters 

The values of the b parameters for the different h(s,b) control structures are tab- 
ulated in table n for 10 -percent increments in control effort. The increase in control 
effort is closely related to increased controller gains. Between the 90 and 100 percent 
of vinity control levels there is a sharp increase in the controller bg and bg gains (for 
the lead-lag and second-order lead-lag controllers). This occurs at the control effort 
levels at which the pressure variation almost vanishes (fig. 4). In the vicinity of this 
rapid change, the most interesting controllers can be chosen. Further analysis will be 
directed to the 90 and 100 percent of unity normalized control effort conditions. 


Effects of Measurement Noise 

The controller designs were obtained with no measurement noise. These controller 
designs were then evaluated for the case of measurement noise to obtain the noise 
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sensitivity of the controllers (fig. 5). 

fi O 

The D/V disturbance to noise ratio of psd values was set at 3. 316x10 [(cm / 
sec-rad)/(N/cm^ -rad )] ^ or 5870 [(in. Vsec-rad)/(lbf/in.^-rad)]^. This resulted in a 
pressure error I ^71^/1 ^ 7 ! v of 0.27 for the 90 percent of unity control effort 
lead-lag controller design, a rather harsh condition. At this noise level the control ef- 
fort was increased 25 percent at the normalization point due to the noise. For compar- 
ison, the control effort values from figure 4 for the case of no measurement noise are 
listed beside each data point in figure 5. With measurement noise the pressure error 
does not drop to vanishingly small values. Bo th the pressure error created by the con- 
troller in attempting to control the noise | P 7 ly and the noise induced control effort 
1 q increase with higher gain controllers. A more detailed design would consider 
different noise levels and optimize the controllers for the noise. 


System Frequency Response 

The frequency response of the system indicates how the cost function averages were 
obtained by giving an "vinaveraged ” view of the controller performance as a function of 
frequency. The frequency response of the P^/d pressure variation to disturbance rela- 
tionship is considered first. The frequency response magnitude | P;y(j(i))/d is plotted in 
figure 6(a) for the 90 percent of vuiity control effort designs. The open-loop response 
(from fig. 3(a)) is plotted again for comparison purposes. The proportional control 
forces the system above the open-loop response at the higher frequencies. The deriv- 
ative control is above the proportional at low frequencies. The lead-lag and second- 
order lead -lag controllers show improvement in performance over the entire frequency 
range. The open- and closed-loop curves must all eventually coincide, but this occurs 
above the frequency range considered here. 

Figure 6(b) is a similar plot of | P,^(ja>)/d| but for the higher gain \mity control ef- 
forts designs. The derivative controller now has more attenuation at the lower frequen- 
cies. The attenuation of all the controllers is generally greater than the figitre 6(a) 
levels. This is consistent with the increased control effort. 

The frequency response magnitude of the suppressor flow to disturbance flow 
|Qa(jo))/d| relationship is considered next. These responses are plotted in figure 7 for 
the controller designs with 90 and 100 percent of \mity control effort. The responses are 
generally flat showing the control effort is about equally effective over the frequency 
range. Variations in magnitude from the nominal level tend to indicate deviations from 
an ideal control. 
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System Eigenvalues 


A study of the system closed-loop eigenvalues, which are the roots of the system 
characteristic equation 1 + G(s) = 0, is made next. The system's two worst eigenvalue 
damping ratios were generally associated with the 2.3- and 24-hertz structural modes. 
The 2. 3 -hertz mode is considered first. The damping ratio for that mode is plotted ver- 
sus control effort in figure 8(a). Dampii^ decreases with decreased control for all de- 
signs, but the derivative curve falls off considerably faster. The 24-hertz mode is con- 
sidered next (fig. 8(b)). Here, derivative action is relatively more effective for the 24- 
hertz mode than it was for the 2. 3 -hertz mode. 

To illustrate the effect on the other modes, all the system structural mode eigen- 
values are plotted in figure 9 for a specific controller, namely, the first-order lead -lag 
90 percent of unity normalized control effort design. With this control, the damping 
ratios are all increased except for that mode associated with the 26. 8 -hertz structural 
mode. The discrete points are connected by lines for illustration only. 


Sensitivity to Structural Frequency 

It is of interest to examine system behavior for changes in the shuttle model. The 
change considered is the variation of the natural frequencies of struct\iral modes. 

The effect on the closed -loop damping for variations of the 24 -hertz structural mode 
is considered first. This mode had about the lowest damping of all modes both with and 
without control. The performance usir^ the 90 percent of imity control designs is con- 
sidered first (fig. 10(a)). The open-loop damping is also plotted for comparison. Judg- 
ing from the cvirves, the choice of a nominal frequency of 24 hertz for the worst -case 
condition was about right. The system is somewhat sensitive to plant changes. The 
more complex controllers do not always achieve the best damping. But all the con- 
trollers maintained system stability. 

The stability of the unity normalized control level designs are considered next in 
figure 10(b). The derivative and second -order lead-lag results were nearly identical for 
this mode. The increase in control effort has not noticeably improved the system sta- 
bility. To show this more clearly, the unity control lead-lag design (from fig. 10(b)) is 
compared with the corresponding 90 percent of xmity design (from fig. 10(a)) in figure 11. 
The lower control level design has about the same average level of stability for this 
mode. 

The sensitivity of the system to changes in the nominal 2. 3 -hertz structural mode is 
shown in figure 12 for controller designs with the 90 percent of unity control design. 
Here, the derivative controller can not maintain system stability for a decrease in nat- 
imal frequency of about 2. 3 hertz; whereas, the proportional and lead -lag designs show 
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good stability and little sensitivity to the plant changes. The vinity control designs (not 
shown) showed little sensitivity to shifts in the nominal 2. 3 -hertz structural mode. 


CONCLUSIONS 

A frequency -domain parameter optimization technique is used in a preliminary de- 
sign study of controllers for active pogo suppression control of the space shuttle. The 
design approach is considered a viable tool for pogo controller design studies. The tech- 
nique is based on minimizing a cost fvmction using a conjugate gradient search procedure. 
The cost function formulated consists of the averaged square variation of the controlled 
variable, the high-pressure oxidizer pump suction pressure P,^, and a weighted aver- 
aged square control. The space shuttle model was a two-pump, six -structural -mode 
linear state variable matrix model. The model represented a worst-case condition that 
was open-loop \mstable. 

The controller design proceeded by designing a series of controllers with increas- 
ingly complex structures. The four controller designs presented were proportional, 
derivative, lead-lag, and second-order lead-lag. The cost fimction averaged square 
values showed the design compromises between controller complexity and level of control 
effort. The supporting design studies included the effects of measurement noise, the 
frequency responses of the system, system eigenvalue damping ratios, and system sen- 
sitivity to changes in structural mode frequencies. 

The system performance with the worst -case model improved for both an increase 
in controller complexity and control effort. Increased control effort was generally re- 
lated to increased controller gains. The two lowest damping ratios increased uniformly 
with increased control. A sensitivity analysis of structural mode natural frequency 
showed that, except in the immediate vicinity of the design point, the damping did not 
necessarily uniformly improve with either increased control effort or increased filter 
complexity. The sensitivity analysis appears useful in the design approach to help select 
a suitable controller. 

The lead -lag controller was a promising structure; its combination of proportional 

2 

plus approximate derivative action achieved good control performance P,^ with reason- 
able control effort. The lead-lag designs required the lowest levels of control to achieve 
system stability. They could attenuate the dist\u*bance error over a broad range of 
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frequencies, and the designs were effective in increasing system damping for both the 
low- and high-frequency structural modes. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, December 16, 1975, 

505-05. 
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APPENDIX - SPACE SHUTTLE POGO MODEL 


This appendix contains a description of the space shuttle oxidizer (liquid oxygen) 
system dynamic model used for this pogo controller design study. The model is bas- 
ically the one described in references 3 and 4. It has been modified to meet the needs 
of this study as well as other studies conducted at Lewis. The model can be divided 
into a structural portion and a propulsion system portion. The structural model is a 
six -mode approximation describing the liquid oxygen (lox) feedline motion. This model 
is based on one developed by the Aerospace Corporation and uses structural model shape 
data presented in references 3 and 4. The propulsion system portion, modified at Lewis 
is a linear, lumped -parameter model of the lox feedlines, pumps, and thrust-chamber 
dynamics. The two portions of the model are joined through consideration of drag and 
momentum force coupling terms. 

The differential equations that describe the combined structural -propulsion system 
model are presented. To obtain the transfer functions required in the pogo controller 
design study, these equations were put in state -variable form. Then the transfer func- 
tions were conveniently obtained using existing computer programs (ref. 5), 


Symbols 

0 

A cross sectional area, cm 

a acoustic velocity, m/sec 

C compliance, cm^/N 

LPOP compliance, cm^/N 

Cb2 compliance, cm^/N 

Cj discharge coefficient, dimensionless 

F force, N 

^ generalized force, N 

f mean mass flow of lox, g/sec 

g gravitational constant, N/g 

i component of vector 

j component of vector 

L inertance, N-sec ycni^ 

I line length, cm 
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1 1 vertical feedline segment 

1‘L horizontal feedline segment 

mj + 1 LPOP dynamic gain, dimensionless 

mg + 1 HPOP dynamic gain, dimensionless 

M generalized mass, g 

2 

P pressure, N/cm 

^ modal pressure coefficient, N* sec /cm*^ 

2 

P,^ HPOP suction pressure variation, N/cm 

2 

p pressure infeedline lump, N/cm 

O 

Q volumetric flow rate, cm /sec 

q generalized displacement, cm 

q volumetric flow rate in feedline lump, cm^/sec 

R resistance, N sec/cm^ 

tb tank base 

X structural displacement along vertical axis, cm 

z structural displacement along horizontal axis, cm 

^ damping ratio 

O 

p lox density, g/cm 

time constant, sec 

o) structural pole natural frequency, rad/sec 

(p modal displacement coefficient, dimensionless 

i// mode shape vector, dimensionless 

Subscripts: 

a actuator 

bt bubble for LPOP 

b2 bubble for HPOP 

c thrust chamber 

d discharge line 

d7 disturbance at station 7 

I inner pump duct 
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i index 
J injector 
j index 
k index 
I lump 

l\ first (vertical) feedline 

12 second (horizontal) feedline 

n structural mode 

pi LPOP 

p2 HPOP 

tb tank base 

Superscripts; 

(x) X direction 
(z) z direction 


System Geometry 

For dynamic modeling pvirposes a simplified lox feedline - pump-chamber geometry, 
developed by the Aerospace Corporation, was used (see fig. 13). The system consists of 
the lox tank connected to a 31. 1 -meter (102 -ft) long, 0. 25 -meter (0. 82-ft) diameter ver- 
tical feedline, an elbow, followed by a 8.47-meter (27. 8-ft) long, 0. 25-meter (0. 82-ft) 
diameter horizontal feedline connected to the LPOP. A short interpump duct connects 
the LPOP with the HPOP, which pumps the lox through a discharge duct into the thrust 
chamber. The suppressor location is at the inlet to the HPOP. 

The Aerospace structimal model defines seven lox system structural elements that 
can move independently of one another. They are the lox tank, the two feedline seg- 
ments, the two elbows, the LPOP, and the interpump-duct - HPOP - discharge -line - 
thrust -chamber combination. As indicated by the velocity coordinates shown in fig- 
ure 13, the tank and the two feedlines can each move in only their axial directions (xj.j^, 
x^ and the elbows and pump assemblies are free to move both horizontally 

and vertically. Propulsion system dynamics are characterized by volumetric flow rates 
Qj and pressures P. throughout the system. The subscripts i correspond to station 
locations; tb, tank base; 1, first elbow; 2, entrance to horizontal feedline; 3, second 
elbow; 4, LPOP inlet; 5, LPOP exit; 7, HPOP inlet; c, thrust chamber. These station 
definitions are the same as those used in references 3 and 4. 
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Propulsion System D5mamic Equations 


The propvilsion system portion is modeled as a lumped parameter system. The ver- 
tical feedline is broken up into nine lumps and the horizontal feedline into three lumps. 
Figure 14 shows a typical lump, or control volume, used in obtaining the feedline equa- 

XI- 

tions. Shown is the i*^" lump, occurring in the horizontal (second) feedline. The ver- 
tical feedline equations, considering conservation of mass and momentum for each lump, 
can be written as 


^ Pi ^t ^Z-^1 • 

qi _ - qi - ; — 






qi = 


^q 


Pi -1 ^Z“^l 


71 


i = 2,3,.. .,9 




Pi = ^^i’Vl^ i = l,2,...,8 


Ps-^i 

Pg = i»2 = -^ (Ql - Q2 + 

Similarly, the following equations can be written for the horizontal feedline; 


k A ^ - PlO , ^2 . ;; 

‘iio ~ S ■ ■ 7“ 'iio ■ 7~ 7“ ^12 


k - 77 ^11 . ^10 , ^1^2 ^ 

qil - - q^j - + + z^2 


^ A A 


R 


I rr 


P12 . Pll ^1^2 • 


qi2 - Q4 - - — qi2 - — + — + — ^12 
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Pio “ ^10 “ 


The cross-sectional areas of the two lines are assumed equal (Aj = Ag = A), thus 
the compliance per lump is 


C 

where g is the gravitational constant, p is the lox density, a is the speed of sound in 
lox, (=pZ/Ag) is the inertance per lump, is the resistance per lump, and is 
the tank base pressure. 

The eqtiation for the pressure at the inlet to the LPOP can be written considering the 
cavitation bubble compliance, the motion of the elbow (3), and LPOP (4) as they affect 
the fluid system. The result is 


_ gA? 

pa^ 


Pl 2 - P 4 — ^ (^4 

’^bl ^bl 


Xs) + 


'bl 


'bl 


where is the cavitation bubble compliance and Ag is the LPOP exit area, which 
equals the interpump duct area. The equations for the interpump duct were written con- 
sidering the duct as one lump with resistance Rj, inertance Lj, and cavitation com- 
pliance ^b2 at the HPOP inlet. The resulting equations are 


Qg - L [(m^ + 1)P^ _ p^ _ (Rj + Rpj)Q5 + + RjAgZ J 

Si 

Pg = (mi + 1)P^ - RpiQg - Lpl^S Si^3^4 


P7 = [Q5 - Qj. + Qa + Q(J7 + (A4 - Ag)z7] 

Si 

where Lpi is the LPOP inertance, mi + 1 is the LPOP gain, Rpi is the LPOP re- 
sistance, is the suppressor flow rate (located at the HPOP inlet), Q^rj is the flow 
disturbance occurring at the HPOP inlet, and A^ is the HPOP suction area. All param- 
eters associated with the pumps and interpump duct were obtained from references 3 
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and 4, Completing the propulsion system model, the equations for the discharge line 
and thrust chamber pressure and flow can be written as 


and 


Q = 


^ + ^p2 


[P? - Pc ■ ^d ^p2^^c \^l\ 



where the following parameters were also obtained from references 3 and 4 : is the 

discharge line inertance, Lj is the injector inertance, Lp 2 is the HPOP inertance, 
is the discharge line resistance, Rj is the injector resistance, Rpg is the HPOP 
resistance, R^ is the thrust chamber resistance, and is the thrust chamber time 
constant. 


Structural Dynamic Equations 

As indicated previously, the structural model used is based on that of references 3 
and 4. In the model seven structural elements are identified that characterize the lox 
system. Primary structural output quantities of interest are the horizontal z and ver- 
tical X velocities of each element. The structural elements and their associated 
indexes are tb, tank base; Zl, vertical feedline; 1, first elbow; 12 , horizontal feedline; 
3, second elbow; 4, LPOP; 7, interpump duct -HPOP -discharge line -thrust chamber 
combination. 

The structural model was developed as was done in references 3 and 4 except for 
the treatment of feedline drag forces. The modal equations can be written as 

+ 2 = -^n " = 1.2 6 

th 

where M„ is the generalized mass of the mode, q„ is the generalized displace- 

ll xy. n 

ment of the n mode, is the n mode damping ratio, is the n mode natural 
frequency, and is the generalized force associated with the n^'^ mode. 

The model in references 3 and 4 has six modes, which are assumed to adequately 
describe the system at a particular flight condition. In this model the structural veloc- 
ities Xj and z. can be related to the six generalized velocities q^^ by 
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n=l 


i = Zl,l,3,4,7,tb 



i = Z2,l,3,4,7 


where and are components of the mode shape vectors. For example, the 

ri_ n 

n mode shape vector is written as 




n 


.j’S'Vn) 

^'( 1 ) 

,pW(tb) 

^^(7) 


where the ordering of structural velocities is arbitrary. The remaining structural out- 
put, tank pressure, can be related to the generalized accelerations, as shown in refer- 
ences 3 and 4, as 
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D 

i=l 


where coefficients ^Ath) are given. 

The generalized force on the n*^" mode was shown in references 3 and 4 to be 




k=l J i i 


where 

i = n, 1,3,4, 7, tb 
j =Z2,1,3,4,7 
n = 1,2, ... ,6 

The first term accounts for forces generated by lox tank out -flow and the second for 
forces due to tank pressure variation. The terms and are the actual forces 
on the structure due to flow rate, pressure, and structural velocities at the seven de- 
fined locations. These forces are made up of (1) drag forces on the feedlines, interpump 
duct, and discharge line, (2) momentum forces at the two feedline corners, the LPOP, 
and the HPOP, and (3) forces on the thrust chamber. 

In a manner similar to that used in references 3 and 4, these forces can be ex- 
pressed as follows in terms of structural velocities and flows and pressures previously 
defined. 

Vertical and horizontal feedlines: The modeling here differs from that used in ref- 
erences 3 and 4 since the feedline models are lumped instead of distributed. The drag 
forces are 
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First elbow. 




= -P2A2 - f (2 ^ - z, 


^ V 


Second elbow: 


LPOP: 


^2 


F("^ = P4A2+f (2 ^ + 23+2x3, 
^2 


Ff = -P4A2-f(2^-HX4j 
^2 


= -P5A3 - ' f Tg " 


Interpump-duct - HPOP -discharge line -thrust chamber: 


- A c P 
^7 ~'^c‘^rc 


F<^> = A3Ri(Q, - AjZ,) + A4R^(Q^ - A^z,) * P,Aj - P^A^ + 2t (_' - ^ U A^P 


i-A 


A3 


Here, f is the mean mass flow of the lox, A^ is the thrust chamber area, and Cj is 
the discharge coefficient. 
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Parameter Values 


The nominal parameters that characterize the propulsion system are given in 
table m. All, except for line inertance and compliance and for actuator frequency, are 
obtained from references 3 and 4. 

Bxibble compliances and C ^2 and pump gains m^ + 1 and nig + 1 are espe- 
cially dependent on flight condition. The values shown in table III define the expected 
parameter ranges. Structural parameters required are mainly the modal vector data, 
plus modal frequencies, masses, and damping ratios. Table IV presents these data, 
which are adapted from data in references 3 and 4, for the two flight conditions of 
interest - end burn and after solid rocket separation. These two conditions were chosen 
because preliminary analysis showed their pogo instabilities to be most pronounced. 


Nominal Worst -Case Model 

A nominal worst -case model was developed to reflect the dynamics of the two most 
critical flight conditions, end burn and after solid rocket separation, and to achieve an 
open-loop model of the more imstable condition. The worst-case model is basically the 
reference 3 model at end burn except for the following: The end burn structural mode at 
2. 81 hertz (El) was replaced by a less stable mode (Al) from the after separation con- 
dition (see table IV). The frequency of the fifth end burn mode was shifted from 26. 8 
hertz to 24 hertz making the system open-loop vmstable. Also, the LPOP and HPOP 
gains were set to the maximum expected values shown in table in, and the pump inlet 
bubble compliances were set to the minimum values shown in table in. 
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TABLE I. - DAMPING RATIOS AND NATURAL FREQUENCIES 
OF OPEN -LOOP PLANT G(s) COMPLEX ROOTS 


Damping ratio 

Natural frequency, 
Hz 

Damping ratio 

Natural frequency, 
Hz 

0.000784 

^2. 31 

0.0190 

^26. 30 

.165 

2. 89 

.0180 

28. 11 

. 00986 

^4. 72 

.0190 

32.99 

. 004 76 

^8. 38 

.0130 

35.91 

.267 

8. 91 

.00750 

39.30 

. 0449 

14. 10 

.00441 

42.26 

.0412 

19. 42 

.00261 

44.63 

. 00898 

^22. 50 

.00162 

46.36 

. 122 

23. 30 

.00110 

47.41 

-. 00151 

^23. 83 




Compare with uncoupled structural mode natural frequencies of 2. 30, 
4.73, 8.49, 22.50, 24.00, and 26. 80 Hz, respectively. 


TABLE n. - CONTROLLER PARAMETERS FOR DIFFERENT 


CONTROL EFFORT DESIGNS 


Control 

effort, 

(^)n. 

percent 

of unity 

Controller structure and parameters 

Proportional 

gain. 

Derivative 

gain, 

^2 

Lead-lag 

Second -order lead -lag 

Gain, 

^3 

Zero, 

b4, 

rad/sec 

Gain, 

Damping 

ratio, 

’^6 

Zero natural 
frequency, 

^7, 

rad/sec 

10 

— 



4.8 

108 

6 

0.43 

237 

20 

(a) 

— 

8.5 

128 

13 

.68 

258 

30 

20.2 

— 

12.0 

122 

22 

. 52 

257 

40 

25.2 . 

— 

14.5 

109 

32 

.48 

251 

50 

30. 2 

— 

17.0 

97 

46 

.44 

245 

60 

35.2 

— 

20.2 

87 

67 

.42 

240 

70 

40.2 

— 

25.0 

75 

100 

.40 

236 

80 

45. 3 

(a) 

35 

62 

160 

.395 

232 

90 

50. 5 

1 

50 

51 

370 

.385 

228 

100 

55. 7 

80 

200 

46 

7000 

.380 

224 


^System unstable due to controller gain being below lower limit for stability. 





TABLE m. - PROPULSION SYSTEM PARAMETERS 


a 

Cbl 

Sb2 

f 

g 

I 

L, 


"pi 

■'p2 


+ 1 

“2 + ^ 
Rr 


p2 


R. 


2 2 

First feedline cross sectional area, cm ; in. 

2 2 

Second feedline cross sectional area, cm ; in. 

2 2 

LPOP exit area, cm ; in. 

2 2 

HPOP suction area, cm , in. 

2 2 

Effective thrust chamber area, cm ; in. 

Acoustic velocity a, m/sec; in. /sec 

Lump compliance, cm^N; in. VlW 

HPOP compliance, cm®/N; in./lbf 

LPOP compliance, cm^/N; in./lbf 

Mean flow of liquid oxygen, N- sec/cm; Ibf- sec/in. 

Gravitational constant, cm/sec ; in. /sec 

Line length, cm; in. 

2 2/5 

Lump inertance, N'sec /cm; Ibf- sec /m. 

2 2 5 

Inner pump duct inertance, N- sec /cm; Ibf - sec^/in. 

LPOP inertance, N-secVc™; Ibf-sec^/in. ® 

2 5 

HPOP inertance, N- sec /cm; Ibf- sec/in. 

2 5 

Discharge line inertance, N- sec /cm; Ibf - sec/in. ^ 

Injector inertance, N- sec /cm; Ibf- sec/in. 

LPOP gain 

HPOP gain 

5 5 

Lump resistance, N- sec/cm ; Ibf* sec/in, 

/ 5 5 

Inner pump duct resistance, N* sec/cm ; Ibf ♦ sec/in. 
HPOP resistance, N- sec/cm^; Ibf- sec/in. ^ 

5 c 

Discharge line resistance, N* sec/cm ; Ibf- sec/in. 

5 / 5 

Injector resistance, N- sec/cm ; Ibf- sec/in. 

/ 5 5 

Thrust chamber resistance, N- sec/cm ; Ibf- sec/in. 

Q 

Liquid oxygen density, g/cm ; Ib/in. 

Thrust chamber time constant, sec 


487; 75.5 
487; 75.5 


201; 31.2 
83.3; 12.6 
1020; 158 
5. 183; 20 400 
5.94; 0.250 


5.94 - 23.8; 0.25-1 
2.38 - 2.98; 0. 1 - 0. 125 
4.01; 2.29 


7.46x10 

2.25x10" 

1.73x10" 

4.33x10" 

2.60x10" 

8.67x10" 


981. 5; 386.4 
- 6 . 


^-6 


345; 136 
1.77X10"^ 
5.36x10"^ 
4. 12x10"^ 
1.03x10"^ 

6. 18x10"^ 
2.06x10"^ 
1.31 - 2.20 
1.42 - 1.54 


-5 


4.14X10""; 9.85x10 
1.91x10"^; 4.53x10"^ 
4. 58x10"^; 0. 109 
8.49x10"^; 2.029x10"^ 
6. 17X10'^; 0. 147 
4. 33x10"^; 0. 103 
1.055; 3.813x10"^ 
4.5X10"^ 




Parameter 


End burn 


After separation 







Mode number, 

1 




1 

2 

3 

4 

5 

6 

1 



Aerospace designation 



El 

E2 

E7 

E30 

E34 

E35 

A1 

Modal displacement 


-0. 156 

-0. 198 



0.222 

0.022 

0. 529 



-0. 330 

-0.245 


-0. 752 

0.370 

0. 0134 

0.208 



-0. 330 

-0.245 

-0.0787 


0. 370 

0. 0134 

0.208 



0.255 

0.239 

-0. 556 

-1.71 

-2. 13 

2. 63 

1.28 



0.255 

0.239 

-0. 556 

-1.71 

-2. 13 

2. 63 

1.28 


(tb) 

-0.209 

-0. 334 

-0.259 

-7. 74 

4.443 

5.336 

0.257 


cp|^'U2) 

0. 0367 

-0. 352 

0.431 

-0.372 

0.017 

-0. 09 

0.094 


crim'd) 

0. 0367 

-0. 352 

0.431 


0.017 

-0. 09 

0.094 



0.0367 

-0. 352 

0.431 

-0. 372 

0.017 

-0.09 

0. 094 


cp|^'(4) 

0. 846 

0. 50 

-1.95 

-0. 12 

0. 99 

-1. 56 

1. 68 



0. 846 

0. 50 

-1. 95 

-0. 12 

0. 99 

-1. 56 

1. 68 

Modal pressure 

^.(tb) 

0 

0 

0 

0 

0 

0 

-2. 16x10"^ 

coefficient, 

XT 2/3 

N- sec /cm 









Generalized mass 

Mi 

292 

292 

292 

292 

292 

292 

1282 

Damping 

^i 

0. 01 

0. 01 

0. 01 

0.01 

0. 01 

0. 01 

0. 01 

Structural pole 

2!ra'^ 

2. 81 

4. 73 

8. 99 

22. 5 

26.8 

27.2 

2. 30 

natural fre- 








quency, Hz 







































3i 







Proportional, 


^ o> a> 

m o J 



N(^) ‘uoi.iBijeA ajnssajd paziieiojON 



N(id) ‘UOUBIJBA ojnssajd pozi|eujJON 







Derivative 



Frequency, Hz Frequency, Hz 

(a) Controller designs for 90 percent normalized control effort. (b) Controller designs for 100 percent normalized control effort. 

Figure?. - Frequency response magnitude of suppressor flow Q, to feedline disturbance flow d. 



Eigenvalue damping ratio 






Eigenvalue damping ratio 



(a) 90 Percent of unity normalized control level. (b) Unity normalized control effort. 

Figure 10. - Sensitivity of eigenvalue damping ratio to variations in 24-hertz structural mode fre- 
quency. 



Figure 11. - Sensitivity comparison of lead-lag controller structure 
with two normalized control levels. 
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Second-order lead-lag 



Figure 12. - Sensitivity ol eigenvalue damping to variation in 
2. 3-hertz structural mode frequency for 90 percent of 
unity normalized control level controllers. 




r- Oxidizer tank 



High pressure oxidizer pump (HPOP)-^ 
Figure 13. - Schematic of liquid oxygen system model. 



Figure 14. - Feedline control volume. (Lump 
length. (; pressure in i'*^ iump (second 
feedline), Pj: flow entering i**^ lump 
(second fe^line). qp horizontal feedline 
velocity, zj 2 i cross sectional area (second 

feedline). A 2 .) 
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